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ABSTRACT. Given a knot K in an Euclidean space R™, and a finite dimensional subspace V C 
C°°(K), we express the expected number of critical points of a random function in V'm terms of 
an integral-geometric invariant of K and V. When V consists of the restrictions to K of homogeneous 
polynomials of degree I on R™, this invariant takes the form of total curvature of a certain immersion 
of K. In particular, when K is the unit circle in R 2 centered at the origin, then the expected number of 
critical points of the restriction to if of a random homogeneous polynomial of degree t is 2V3£ — 2, 
and the expected number of critical points on A" of a random trigonometric polynomial of degree k is 
approximately 1.549/c. 



To the memory of my mathematical hero, Vladimir Igorevich Arnold 

Contents 



Introduction 1 

1 . An abstract result 2 

2. Polynomial Morse functions on knots 6 

3. An application 8 

4. Critical points of trigonometric polynomials 11 
References 12 



Introduction 

A celebrated result of Fary and Milnor [2, 5] states that the expected number of critical points 
of the restriction to a knot K <— > M 3 of a random linear map h : R 3 — > E is equal to an integral- 
geometric invariant of the knot, namely, its (suitably normalized) total curvature. It is natural to 
ask how this result changes if, instead of random linear maps, we look at random homogeneous 
polynomials P : M 3 — > M of a given degree I. It is convenient to investigate an even more general 
situation. 

Suppose that U is an oriented real Euclidean space of dimension n and K U is a knot in U, 
i.e., an smoothly embedded S . We assume that K. For any positive integer I we denote by 
Qe(U) the space of symmetric ^-linear forms on U. Each such form defines a polynomial function 

P^-.U^R, P$(u) = $(«, ...,u). 

i 

We denote by /$ the restriction of P$ to the knot K. We will refer to such functions as polynomial 
functions on K of degree i. For random the function /$ is Morse, and we denote by ^ K {^) its 
number of critical points. We can regard ^k(&) as a random variable and ask how is its expectation 
related to the global geometry of K. 
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To describe this relationship, denote by S the unit sphere in Qe(U) with respect to the natural 
metric (—,—}e induced by the Euclidean metric on U. The expected number of critical points of 
is the real number f/ K defined by 



area (b l ) Jge 



In Theorem 2. 1 we express \x l K as an integral geometric invariant of K. This invariant is described in 
terms of the Veronese map Vt : U — > Qg(U) uniquely determined by the equality 

(<s>,V e (u)) e = $(u,...,u), Vuet/, $ e Qe(U). 

i 

The restriction of the Veronese map to K produces an immersion 

V e : K -> Q £ (U) 

that we called the t-th Veronese immersion of K. In Theorem 2.1 we show that f/ K is the total 
curvature of the ^-th Veronese immersion. The case t = 1 of this theorem is the celebrated result of 
Fary and Milnor [2, 5]. 

As an application we compute i/ K , when K is the unit circle in M. 2 centered at the origin. More 
precisely, we show that in this case \i l K = 2\/3£ — 2. In other words, the expected number of critical 
points of a polynomial function on K of degree I is 2^J?>i — 2. 

We obtain Theorem 2.1 as a special case of Theorem 1.1 that describes the expected number of 
critical points of a random function belonging to a fixed finite dimensional subspace V C C°° {K ) 
satisfying a nondegeneracy condition (1.1). 

Theorem 1 . 1 has another interesting consequence. More precisely, in Theorem 4. 1 we show that 
the expected number of critical points on the unit circle of a random trigonometric polynomial of 
degree n is 

_ 5/3 3 (n + l)' 

where (3k denotes the k-th Bernoulli polynomial. In particular 

^3 



2> 



T n ~ 2ny - as n — > oo. 

Notations. We will denote by cr n the "area" of the round n-dimensional sphere S n of radius 1, and 
by u n the "volume" of the unit ball in IR n . These quantities are uniquely determined by the equalities 
(see [7, Ex. 9.1.11]) 

r(i/2) n 



r„_i = nu n , ijj n = — — , r(l/2) = V7T, 

' r(l + n/2)' v 1 ' 



where V is Euler's Gamma function. 



1. AN ABSTRACT RESULT 

Let U and V be a real, oriented Euclidean space of dimension n. We denote by (— , —)u the inner 
product in U and by S(U) the unit sphere in U. Suppose that K ^ U is a smooth knot in U, i.e., 
the image of a smooth embedding S 1 U. Fix an arclength parametrization of K 

s i-)- x(s), < s < L, 

where L denotes the length of K. 
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Assume that we are given a finite dimensional subspace V C C°°(U) of dimension N satisfying 
the nondegeneracy condition 

VsG[0,L], 3veV: d x{s) v(x' (s)) + 0, (1.1) 

where d x v denotes the differential of the function v at x. We fix an Euclidean inner product (— , — )v 
on V and we denote by S(V) the unit sphere in V centered at the origin. 

As explained in [7, §1.2], the condition (1.1) implies that for generic v G V, the restriction of the 
function v to K is a Morse function. We denote by hk{v) its number of critical points. Note that for 
any nonzero scalar t G M the restriction of tv to K is Morse if and only if the restriction of v to K 
is such. Thus, it suffices to concentrate on functions v G S(V). The main goal of this section is the 
computation of the expected number of critical points of a random function in S(V), i.e., the quantity 

^■■= T^TVT / VK(v)\dS(v)\ = \* N ^{S(V)) [ fi K (v)\dS(v)\, 

area {S(V) ) J S (V) I Js(V) 

where \dS\ is the "area" density on S(V). 

To describe the result we need to introduce the main characters. Set V y := Hom(V,M) and 
observe that we have a smooth map £ : K — > V y that associates to each point x(s) G K the linear 
map 

& : V -> R, v ^ d x[s) v(x'(s)) G R. (1.2) 
Using the metric induced isomorphism V y — > V we obtain a dual map 

: K -> V. (1.3) 
The nondegeneracy condition (1.1) implies that ^ 0, Vs. We can thus define a smooth map 

u:K^S{V), := f4&- 

ls«l 

Theorem 1.1. Le? J7, K and V be as above then 

lil ■= — / ^(«) |rfS(«)| = - / KWI \ds\. 

CTN-l JS{V) n JK 



Proof. Set 

E v := {(x,v) G K x 5(V); i/(as)±w}. 
Note that can be alternatively defined as the zero set of the function 

F : K x 5(F) -> E, F(sb, «) = (i/(a?), w) V) 

and the differential of F is nonzero along the level set {F = 0}. This shows that E v is a smooth 
submanifold of K x S(V) and 

dimF^ = dim S(V) = N - 1. 

We denote by <?£ the metric on Fj, induced by the natural metric of K x S(V), and by |cZVe| the 
associated volume density. We have two natural (left and right) smooth maps 

S{V). 

The fiber of A over x G F is the Equator 

Fa, := {i> G S(V); v ± u(x) }, 
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while the fibers of p are generically finite. Tautologically, the restriction of A to any fiber of p is 
injective. More importantly, for any v £ S(V) the subset A(p _1 (i>)) C K is the set of critical points 
of the restriction of v to K. Hence, for generic v G V we have 

Hk(v) = #p~ l {v). 

We have thus reduced the problem to computing the average number of points in the fibers of p in 
terms of integral-geometric invariants of the map v. We will achieve this in (1.10). 
The area formula (see [3, §3.2] or [4, §5.1]) implies that 

#p- 1 (w)\dS(v)\= [ J p (x,v)\dV E (x,v)\, (1.4) 

S(V) JE 

where the nonnegative function J p is the Jacobian of p defind by the equality 

P *\dS\ = J p -\dV E \. 

To compute the integral in the right 4iand side of (1.4) we need a more explicit description of the 
geometry of E v . 

Let (xq, vq) G E v . Fix an orthonormal frame e\, . . . , of the tangent space of E XQ at vq. 
Assume xq = x(sq). Then the tangent space of E u at (xq, Vq) consists of tangent vectors 

sx'(so) ®v e T X0 K ®T V0 S(V), s£R, 

such that 

s{u'(s ), v )v + (v(s q ),v)y = 0. (1.5) 

Define 

e = x'(s) - (v'(s ),v ) v u(s ). 

For simplicity we set 

H(x ,v ) := | (u'(s ),v )v |- 
The equality (1.5) implies that the collection 

eo, ei, . . . , 6 TV— 2 

is an orthogonal basis of T^ X0 ^E V . Moreover, the length of eo is 

|e | = ^1 + p{x ,v ) 2 . 
If we denote by dS x the area form on E x then we see that at (xq,vq) we have 

ds A dS x (e , ei, . . . , eAr-2) = ±1. 

Hence at (xq, Vq) we have 



\dVE\ = I e 1 • \ds A d5a;| = \Jl + p(x ,v ) 2 \ds A d5a;|. (1.6) 
The differential of p at (xq,vq) is the linear map 

P* '■ T(x ,v )Ev -> T V0 S(V) 

given by 

eo ^ (v'(s ),v )v ■ ^(so), ei 1-^ e i5 1 < i < iV - 2. 

We conclude that 
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Using (1.6), (1.7) and the co-area formula for the map A [6, Prop. 9.1.8] we deduce 



J p (x,v)\dV E (x,v)\= / fi(x,v)\dS x (v)\ )\ds(x)\. (1.8) 

E JK \J E x J 

v v ' 

=:J(x) 

To proceed further we need the following elementary result. 

Lemma 1.2. Suppose W is an (m+l)-dimensional oriented real Euclidean space with inner product 
( — , — ), and Vq € W. Denote by S(W) the unit sphere in E, and by \dS\ the "area" density on 
S(W). Then 



I ( r o) = / \(v ,w)\\dS{w)\ = u m \v \. 
JS{W) 

Proof. Fix an orthonormal basis (eo, e±, . . . , e m ) of W and denote by (wo, ■ ■ ■ , w m ) the resulting 
coordinates. Observe that for any orthogonal transformation T : W — > W we have I(Tvq) = I(vo) 
so that I(vq) depends only on the length of i>o. Thus, after an orthogonal transformation, we can 
assume v$ = ceo. We can then write 

I(ce ) = \c\[ \(e ,w)\\dS(w)\ = \c\ [ \w \\dS(w)\. 
JS(W) JS{W) 

Denote by S + (W) the hemisphere wq > 0. We have 



/(ceo) = 2|c| / wo\dS(w) 
Js+(w) 



The upper hemisphere S+(W) is the graph of the map 



w = yl - r 2 , r := \J w 2 -\ h w 2 n < 1. 

Observe that 

r 2 1 1 
\Vw \ = . and \ dS ( w )\ = , * Wm\ = —\dV m \, 

where | dV m \ denotes the Euclidean volume density on the m-dimensional Euclidean space with co- 
ordinates (wi, ■ ■ ■ , w m ). We deduce 



I{ce ) = 2\c\ / \dV m \ = 2\c\u r 

Jr<l 



□ 



Using Lemma 1.2 in the special case when m + 1 = N — 1, W is the hyperplane containing the 
Equator E XQ , x$ = x(sq) € K, and vq = u'(so) we deduce 

J(x ) = I(u'(s )) = 2|l/'(s )|<^Af-2- 

From ( 1 .4) we now deduce 

0~N-1 Js(V) 

Using the equalities 

r(i/2) r 



/ #p-\v)\dS(v)\ = / \v\s)\\ds\. (1.9) 

JS(V) JK 



<t N -i = Noj n , u r 



T(l + m/2)" 
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we deduce 

2cjjy„ 2 _ 2 T(l + N/2) _ 1 
o-jv-i ~ JVT(l/2) 2 r(JV/2) ~ vr' 
where at the last step we have use the classical identities 

r(l + N/2) = jT(N/2), r(i/2) = v^F- 



We have thus proved 



/ #p~ 1 (i,)|d5'(v)| = - / \v'(s)\\ds\. (1.10) 



°"iV-l J5(V^ 



2. Polynomial Morse functions on knots 



□ 



Let U and K U as above. We want to apply the results in the previous section to a special 
choice of V. Fix a positive integer I. Denote by Qe(U) the space of symmetric ^-linear forms on U, 
or equivalently, the space of homogeneous polynomials on U of degree L This will be our choice of 
subspace V C C°°(U). The nondegeneracy assumption (1.1) translates into the condition K. 
Moreover, 

dimQ,(C/)= ( n + £ £ " r 

The Euclidean metric on £/ induces an inner product (— ,— ^onQ^C/). Denote by §> e the unit sphere 
in Q^(CJ). For any <I> G we obtain a function 

/* :iT^R, /*(aj) = $(^_^). 

The function /$ is a Morse function for almost all $ G S £ . We denote by ^k(^) the number of 
critical points of fp. The main goal of this section is to describe the average 



r N e -l J$e 

in terms of integral-geometric invariants of K. We will achieve this by reducing the problem to the 
situation investigated in the previous section. 

We begin with a linear algebra digression. The £-th Veronese map is the linear map 

V = V e : U®- ■ ■®U ^ Q e (U), 

s v y 

I 

uniquely determined by the requirement 

($,V(u 1 ,...,'Ui)) e = *(«!,..., u e ),, Vu 1 ,...,u e eU, $eQ*(l7). (2.1) 
Note that for any smooth paths 

t h> Ui(t) G U, i = l, ...,£, 

we have 

— V(m, ...,u e ) = V(u!,u 2 , ... ,u t ) H h V(«i, . . . ,u e _i,u e ), 

where a dot indicates i-differentiation. For any u G U we set 

M u := V(tt,...,tt). 
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With these notations we deduce that 

fc(x) = {$,M x )t, \/x G K, $eS £ . 
We deduce that xq = x(sq) is a critical point of /$ if 

where a prime ' indicates ^-differentiation. Observe that 

M' x{s) =£V(x'(s),x(s),...,x(s)), 

i-l 

Since x(s),x'(s) ^ 0, Vs we deduce 

K{s) + °> Vs - ( 2 - 2 ) 

Define 

,/(,) = -±-M' a(a) . 

I !B(«)I 

Using (1.10) we deduce 

A = -\ W^)\\ds\. (2.3) 
To formulate this in a more geometric fashion observe that (2.2) implies that the map 

K 3 x ^ Vix^^x) G Q £ (2.4) 

is an immersion. We will refer to it as the i-th Veronese immersion of K. For every s the unit vector 
v t{s) is the unit tangent vector field along this Veronese embedding. The integral in the right-hand 
side of (2.3) is then precisely the total curvature of the ^-th Veronese immersion of K. We denote it 
by Tg(K). We have thus proved the following result. 

Theorem 2.1. Let K CU be a knot in the Euclidean space U such that K. Then for any positive 
integer I we have 

/4r = ~n(K), 

7T 

where Te(K) denotes the total curvature of the i-th Veronese immersion of K defined by (2.4). □ 

Remark 2.2. (a)The case £ = 1 in the above theorem is precisely the celebrated result of Fary and 
Milnor, [2, 5]. 

(b) The arguments in the proof of Theorem 2.1 show something more. More precisely, if X C U 
is a compact submanifold of U, then the expected number of critical points of P$\x is equal to the 
expected number of critical points of the function A o Vg\x, where A : Qe — > R is a random linear 
function of norm 1. The arguments in Chern and Lashof, [1], show that this number is the total 
curvature of the £-th Veronese immersion : X — > Q^(C7). □ 
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3. An application 

We want to compute the higher total curvatures rg(K) when K is the unit circle in U = R 2 
centered at the origin. Denote by Ki the £-th Veronese immersion of K. We begin by providing a 
more explicit description of the Veronese map 

V e : U ® ■ ■ ■ ® U -> Q e (U) 

t 

For any ^-linear form 

* : U x • • • x U -> R 
we define its symmetrization ) € Qe{U) by 

(*>(«!,...,«/) := - ^ VPfu^i),..., 

where denotes the group of permutations of {1, . . . , £}. If u\, . . . , U£ £ U,we denote by u\ the 
linear functionals 

u\ : U ->-R, it- (it) = (ui,u)u. 

Then 

V^(wi, . . . = ( w{ <8> • • • <8> u\ ). 

Denote by (e±, e%) the canonical orthonormal basis of 1R 2 and by (e , e 2 ) its dual basis. An orthonor- 
mal basis of is given by the monomials 

= Q ( {e l p ® (e 2 )^(^) y < i < £ 



For any a; £ U we have 



where 



V / (x) = J2^( aJ )^ > 



j=0 

1/2 



^•(x) = (V/(X), $j)t = Qjix, . . . , x) = Q xja^'. 

The map has an important invariance property. Observe that we have a linear right action of 
SO(U) on Q e 

Q e {U) x 50(17) 9 (S,fl) i-> £ Q*(17), 
where : Q^(J7) — >■ Qe(U) is the isometry defined by 

...,«/) = . . . , Ru t ), 

for all <3? £ and u%, . . . , ui £ J7. The equality (2.1) implies that for any «x, . . . ,ug £ 17 we have 

V e (Ru 1 ,...,Ru e )=7 j R (V e (u 1 ,...,u e )), (3.1) 

where : Qe(U) — > Q^(17) denotes the adjoint of T#. 

Consider the standard parametrization 6 ^ x{9) = (cos 6, sin 0) of the unit circle in U. The £-th 
Veronese immersion of the unit circle is given by the map 

V : [0,2tt] B 9 ^ ( V o (0), V e (6) ) £ R £+1 

/A 1/2 

V$(0) = ( J (cos#) J (sinfl)^, < i < I. 
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We have 



and 



V'{6) = £V e (x'(9), x(6), x(8) ) , , x"(8) = -x(8). 



V"(9) = £V e (x"(9), x(9), x(9) )+£(£- l)V e (x'(9),x'(9),x(9), x(9)). 
If we denote by Rg the counterclockwise rotation of R 2 of angle 9, then 

x{9) = Rox(0), x'{9) = R e x'(0). 

Using (3.1) we deduce that 

V'(9)=tf R V'(0), V"(9) = ^ R V"(9). 
In particular, since the map 7 R : Qe(U) — > Qe(U) is an isometry we deduce 

\v'{9)\ = \v'(o)\, \v"(9)\ = \v"(o)\, ye. 

Observe that for < j < £ we have 

1/2 

(-j (cos 9) (sin 9) ( e - j+ V + {£- j) (cos #) (sin 9) 
j (cos 9)( j ~V (sin 9)^i~V ( -j sin 2 9 + {£ - j) cos 2 9 ) 



V'{8) 



1/2 



(cos 9) (sin 9) ( e - j ~ 1] (- j + £ cos 2 



Next we observe that 

V}(0) = -£(cos9) e - 1 sm9, Vg(0) = £(sm9) £ ~ 1 cos9. 

This shows that for 9 = only the component V^ / _ 1 (0) is non trivial and it is equal to £ 1 / 2 . This 
proves that 

\V'{8)\ = £ 1 ' 2 . 

Thus the arclength parameter along the curve 9 \-t V{9) is s = £ l / 2 9, the length of Kg is 2-k^I 2 and 

± = f-l/2± 

ds dB' 

Next we compute \V"{0)\. If £ - j - 1 > 2, i.e. j < £ - 3 we have Vj'(O) = 0. 

^'(o) = -/, d(o) = o, y;: 2 = 2Q 1/2 . 



Hence 



Thus 



= £ 2 + i 



1 



2 



3^ 2 - 2£ 



rt(K) = - / 
71" Jo 



d 2 



ds = 2V3£ - 2. 



3^-2 



We have thus proved the following result. 
Theorem 3.1. If K is the unit circle in M 2 centered at the origin, then 

f/ K = 2VS£^2, W > 1. 



□ 
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Remark 3.2. (a) The case i = 1 of the above theorem is obvious. For the case £ = 2, i.e., the equality 
H 2 K = 4, we can give two alternate proofs. For the first proof we observe that 

V 2 (cos#,sin#) = ( x9), y(9), z{6) ) (sin 2 9,^2 sinfl cos 9, cos 2 9). 

Observe that the image of the 2nd Veronese immersion is the intersection of the sphere x 2 + y 2 + z 2 = 
1 with the plane x + z = 1 which is a circle of radius Moreover, the Veronese immersion double 

covers this circle. These facts imply easily that [± 2 K = 4. 

For the second proof we observe that the typical pencil of plane conies 

C t = {ax 2 + 2bxy + c 2 y 2 = t} 

has four points of tangency with the unit circle. 




FIGURE 1 . The second Veronese embedding of the circle of radius 1 with center (120, 0). 

(b) The invariant jjir K is invariant under rotations, but not under translations. For example if K is 
the circle of radius 1 in the plane centered at the point (3,0). Then the second Veronese immersion 
of K is given by 

[0,2vr] 3 m ((3 + cost) 2 ,\/2(3 + cost)sint,sin 2 t) G M 3 . 
MAPLE aided computation shows that its pointwise curvature is 



k(S) 



108 (cos {t)y - 54 (cos (t)Y - 360 cos (t) - 238 
-2 (9 (cos (t)f - 19 - 6 cos (£) N 3 



while its total curvature is 

H 2 K 2.065 < 2V3 -2-2 = 4. 
Numerical experimentation suggest that if K{d) is the circle of radius 1 centered at the point (d, 0) 
then 

Jim ^=2. 

In view of Fenchel theorem, this would indicate that for d large the second Veronese embedding of 
K(d) is very close to a planar convex curve. The MAPLE generated Figure 1 seems to confirm this, 
since the second Veronese immersion is contained in a box of dimensions 1 x a x b, where a > 300 
and b > 400. □ 
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4. Critical points of trigonometric polynomials 

As another application of Theorem 1 . 1 , we want to compute the expected number of critical points 
on the unit circle of a trigonometric polynomial of degree < n. 

Let K denote the unit circle in R 2 centered at the origin and denote by V n the vector space of 
trigonometric polynomials of degree < n on K, i.e., functions / : K — > R of the form 

n 

f(9) = ao + ( a^. cos kO + bk sin kO ) , 
k=i 

equipped with the inner product 

(f,9)v = ^- r f(6)g(0)dB. 
27r Jo 

k=l 

eo = 1, efc = V^coskO, f k = V2smk9, k = l,...,n. 

The function 

£:K^V n , [0,2tt] 9 
defined by (1.2) +(1.3) is given in this case by the map 

n 

s \-t £l = k{~ ( sm ks)ek + (cos ks)f k ). 

fc=i 

Observe that 



Thus 

n 

ao + ( a k cos k-9 + bk sin kO ) 

k=l 

We obtain an orthonormal basis of V n 



2 
v fe=l 



We set 



Then 



so that 



fc=l 

. / n \ -V 2 n 

^(s):=-r4 t= ^ fc2 J2k(-(smks)e k + (cos ks)f k 

\^ s \ \fc=i / k=i 

^ ( s ) = (X] fc2 ) ^ fe2 ( ( cosks ) e k ~ (sinks) f k ), 

\k=l J k=l 

/ n \ - 1 / 2 / n \ 1 / 2 

w'(s)\ = (J2 k ) -IE* 4 ) 



Using the classical identity 



\fc=i / \fc=i 



n 1 



k=l 

where is the £-th Bernoulli polynomial we deduce 
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where 



X 



-x 2 + -x, /3 5 (x) 



2 3 



-x. 



2 2 

The equality (4. 1 ) and Theorem 1 . 1 imply the following result. 

Theorem 4.1. The expected number of critical points on the unit circle of a trigonometric polynomial 
of degree < n is 



' 3/3 5 (n + l) 
5/3 3 (n + l) 



In particular 



T n ~ 2n\/ — as n — > oo. 



□ 



Remark 4.2. The above theorem indicates that for large k, a trigonometric polynomial of degree k is 

expected to have approximately 2y|A: 
that 2k critical points on the unit circle. 



expected to have approximately 2 J |& pa 1.549 A; critical points. In particular, it will likely have less 




Figure 2. A "random" trigonometric polynomial of degree 4. 

In the MAPLE generated Figure 2 we depicted the graph of a "random" trigonometric polynomial 
of degree 4 

1.2 cos t + 2.35 sin i — 3.17 cos 2t + 2.71 sin 2t + 1.53 cos 3t - 4.17 sin 3t 

+ cos4t - 1.15 sin At. 

We notice that it has 6 critical points, which is very close to the expected value T4 6.87. □ 
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